Anomalous magnetic splitting of the Kondo resonance 
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The splitting of the Kondo resonance in the density of states of an Anderson impurity in finite 
magnetic field is calculated from the exact Bethe-ansatz solution. The result gives an estimate of 
the electron spectral function for nonzero magnetic field and Kondo temperature, with consequences 
for transport experiments on quantum dots in the Kondo regime. The strong correlations of the 
Kondo ground state cause a significant low-temperature reduction of the peak splitting. Explicit 
formulae are found for the shift and broadening of the Kondo peaks. A likely cause of the problems 
of large- TV approaches to spin-i impurities at finite magnetic field is suggested. 
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The physics of a magnetic impurity in a metallic host 
has been studied for many years, and the nature of the 
ground state and lowest excited states are now under- 
stood through the numerical renormalization group |^] 
(NRG) and the Fermi-liquid theory of the Kondo fixed 
point. The theoretical situation is much less clear for 
excitations away from the Fermi surface; however, it is 
precisely these excitations which are probed by two kinds 
of current experiments. The nonequilibrium conductance 
through a quantum dot in a magnetic field H , believed 
to be the clearest sign of Kondo physics in such systems, 
is determined by the electron spectral function on the 
dot away from the Fermi energy. The Kondo peak in the 
spectral function at the Fermi level splits for iJ > into 
two peaks which move away from the Fermi level with in- 
creasing H . The subject of this paper is an approximate 
calculation of this spectral function, based on the exact 
Bethe-ansatz solution 1^,^ of the s-d model. 

The differential conductance through a quantum dot 
at finite bias but H = was calculated by Meir, 
Wingreen, and Lee using the non-crossing approxi- 
mation (NCA) [0 applied to the Anderson Hamiltonian 

H ^k,acl ,^Ck^a + ^dcli aCd,a + 

UUd.^^rid.i + ^(Vfe,cr4,crC<i,cr + h.C.) (1) 

The NCA is a large- technique which sums a class of 
diagrams sufficiently general to resolve the Kondo tem- 
perature Tq. However, at finite magnetic field the NCA 
was found to give spurious peaks in the density of states, 
and thus a perturbative equations-of-motion (EOM) ap- 
proach was used in which the energy scale Tq does not 
exist (i.e., is effectively zero). The resulting splitting of 
the Kondo-like peaks is exactly twice the s = ^ Zee- 
man splitting. One motivation for this work was to de- 
termine whether the strong correlations of the Kondo 
ground state cause corrections to the peak splitting. 

We develop a nonperturbative estimate of the Kondo 
peak splitting and shape for all values of the ratio 
g^H/To. The decreased splitting can be pictured as re- 



sulting from the oscillation of the impurity spin due to 
its interaction with the lead electrons (the time scale for 
this oscillation is Tq~^ for H — 0). 

We use the Bethe-ansatz solution to calculate exactly 
the density of states for "spinon" excitations, which are 
the only excitations at energies up to g^H from the 
Fermi level. 1^ The density of states for spinon excita- 
tions gives an estimate for the electron impurity spectral 
function. This approximate technique is known at _ff = 
to give results consistent with Fermi-liquid perturbation 
theory and the NRG. ^ Analytic expressions are derived 
for zero temperature, and we comment on the expected 
changes when T > 0; the results on spinon excitations 
given below can in principle be combined with existing 
numerical methods for the finite-temperature thermody- 
namic Bethe ansatz. [l0| ] There are simple limiting forms 
for the peak location and broadening when gfiH 3> Tq. 
In this limit the EOM approach of ||] slowly becomes 
valid as the Kondo energy scale disappears. 

The Bethe-ansatz solutions of the s-d and Anderson 
models give exact results for thermodynamic properties 
such as magnetization and specific heat. Dynamical 
quantities such as the electron Green's function on the 
impurity site in the Anderson model are not known ex- 
actly from the Bethe ansatz (or from any other technique) 
but are required to understand tunneling and other ex- 
periments. For example, the Landauer-type formula for 
nonequilibrium conductance [ pj] depends on the spectral 
function pd,a{^) — (~l/''')IinG'J^(w) over all energies 
between the source and drain voltages. 

The spectral function can be written as a sum over 
energy eigenstates 

Pd,a{uj) = Y [mCd,aW)?5{uJ - {En - Eo)) + 
n 

\{0\4jn)\Hic. + iEn^Eo))] (2) 

At zero magnetic field, an accurate description of the 
shape of the Kondo peak in p{lu) follows from assuming 
the matrix elements to be some constant C for all states 
connected to the ground state by a single spinon (defined 
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below), and for other states. This is successful for two 
reasons: because the Kondo peak in the spectral function 
arises from the Kondo peak in the DOS rather than from 
structure in the matrix elements, and because spinons are 
the excitations onto which the electron operator projects 
most strongly. At H > there is the additional advan- 
tage that all excitations other than spinons acquire a gap 
g^iH and do not affect the spectral function below this 
energy. The value of the constant C, which determines 
the overall weight in the Kondo peak, can be found for 
the symmetric model (e^ = —U/2) from the value of the 
spectral function at the Fermi level pd.AO) = I/ttA, with 
A = Trpo X] l^fcP bare linewidth. |]12| 

The DOS in the exact solution of the Anderson model 
consists of a Kondo peak near zero energy from spinons, 
and a peak at energy U from "holons" (the simplest ex- 
citations in the charge sector). In the Kondo regime, the 
Anderson model reduces to the s-d model (an impurity 
spin with interaction term Jsitp'' aiip, tp the conduction 
electron field) if we are concerned with energies much 
less than U above the ground state. Hence, to under- 
stand the splitting of the Kondo peak for < gfJ,H ^ U, 
the spin excitations are the relevant ones and the s-d 
model is appropriate. Mixed- valence behavior requiring 
a full Anderson model description is expected to begin 
at gfiH > ([/A)^/^, which is a very high field for current 
experiments on quantum dots. 

The exact solution of the s-d model contains a trivial 
charge sector and a spin sector with eigenstates deter- 
mined in the thermodynamic limit from solutions of lin- 
ear integral equations. Just as in Fermi liquid theory, 
a state is labeled by occupation numbers n(i) = 0, 1 on 
the "rapidities" (here the ordered rapidities Xi < Ai+i 
play the role of momenta in the Fermi liquid theory). 
However one complication from ID interaction is that the 
position of rapidities depends on the occupation numbers 
n(i). In the thermodynamic limit, we can introduce the 
density of rapidities Dr{X) — Ai/{Xi+Ai — Xi) to describe 
their distribution. Now a state can be labeled in terms 
of a function cr(A) — n[i{X)]Dr{X) on the real line (where 
i{X) is the integer such that A^ is closest to A). Thus cr(A) 
can be interpreted as the density of occupied rapidities 
A in the state. The total energy and magnetization are 
expressed as weighted integrals of a. All rapidities are 
filled in the ground state, i.e. n{i) = 1 or cr(A) = Dr{X). 

The Bethe-ansatz solution is universal in the limit 
bandwidth L» ^ oo, coupling c = 2J/(1 - 3j2/4) 
with To = £'exp(— 7r/c) finite. The Kondo tempera- 
ture Tq is defined as the halfwidth of the DOS peak at 
T = 0, and is related to Wilson's Tk by Tq = Tk/w, 
w = exp(c-|- l/4)/7r^/2 = 0.41071 . . ., c Euler's constant. 

At zero magnetic field, the energy to add a "hole" 
at Xh {i.e., to set n[i{Xh)] = 0) is tan~-^(exp(7rAft/c)), 
and such a hole has spin i. In a magnetic field, fill- 
ing the lowest rapidities with holes lowers energy, i.e. 
n{i) = for rapidities Xi in the interval {—oo,B), with 



exp(7rB/c) = y/e/2^{g^iBH /Tq). With B > -oo, there 
are two types of spinons: hole-like excitations at A > i? 
with AS — J and particle-like excitations at A < i? 
(where n[i{X)] is changed from to 1) with AS" = — ^. 

The density of states for spinon excitations is Ds{e) = 
Dr{X){dE/dX)~^ , evaluated at the value of A with 
E{X) = e. The density of rapidities Dr for the ground 
state (note that Dr depends on the occupation numbers 
n{i)) was found for nonzero H in the calculation of the 
impurity magnetization, so all we need is the energy of 
the spinon excitation at rapidity A. The DOS Ds{e) con- 
tains both the impurity contribution and a (constant) 
conduction electron part, ignored henceforth. The spinon 
consists of a single i5-function change in the density a (due 
to the change in the occupation numbers n{i)) plus the 
(regular) backfiow (due to the change in Dr caused by 
the change in n{i)). The backfiow (t'{x) for a (5-function 
introduced at Ah > i? satisfies 

/•oo 

a'(A)+ / X(A- A2)a'(A2)dA2 = X(A- Ah) (3) 



with K{x) ~ c/tt{c^ + x^). With B ~ — oo this becomes, 
in Fourier space, d'{p) — 1/(1 + 6^1^'), so the total density 
change is 5(7{p) — a'{p) — 1 = —1/(1 + exp(— c|p|)). p^ ] 
With a magnetic field, B > —oo and (||) can be solved 
by the Wiener-Hopf method. The positive backfiow cr'(A) 
now includes both particles (A > B) and holes {X < B). 
The energy and magnetization of the spinon can be cal- 
culated equivalently by summing over particle or hole 
states. The original hole is dressed by other holes, in- 
creasing its interaction energy but also its magnetization. 
The total change in energy as Ah — > B^ is exactly suffi- 
cient to render the excitation gapless. 

The backfiow for a hole removed at Ap < S satisfies 



a'(A) 



/ K{X- X2)a'{X2)dX2=-K{X- X,). 
Jb 



(4) 



The resulting "particle-like" excitation is gapless as Ac 
B^ and connects smoothly to the hole excitation. For 
Ac ^ B, (j'{X) = —K{X — Ac) and the total number of 
holes removed in the backfiow is 1. 

The Wiener-Hopf technique gives explicit forms for the 
Fourier transforms of 0±(A) = a{X -t- B)9{X + B). The 
results involve the kernel factorization 

^-(^)-(^-(-^»;^-f(bW^ 

ZTT 

exp[-i.-E(l + — log(-a; + i(5))]. (5) 
For a hole excitation (Ah > B): 

""^^'^^^L 2.K^icp/2.)[t-^p) ■ 

For a particle- like excitation (Ae < B): 
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<^_(0) = 1- 



0_(i7r/c) 



dt 



sin(ct) 
7ri^_(0)t' 



^=)*ir_(ict/2^), (7) 



sin(c<)ii:_(ict/27r)e~('^-^=)* 



7rK^{i/2){t~7:/c) 



The total energy for either type of excitation is given by 
E{X) = roe"^/^(e"^/^ + 0-(i7r/c)) 

.,,i7(y^(e-/^ + 0_(^))-i±|^). (9) 

The ground-state density of rapidities crgs(A) also sat- 
isfies an equation of Wiener-Hopf type with source term 
2c/7r(4A^ -f- c^). The solution for Cgs is always peaked 
near A = 0, with crgs(A) = l/2ccosh(7rA/c) for S/c <C 
and crgs(A) = 2c/tt{A\'^ + c^) for B/c > 0. The explicit 
solution combined with (j^) gives the single-spinon 
density of states. The full (thermodynamic) density of 
states up to energy g^H is a sum of convolutions of the 
single-spinon result. The single-spinon density of states 
at zero energy gives the exact spin susceptibility and spe- 
cific heat, whose ratio (the Wilson ratio) does not change 
with magnetic field, p^ ] 





FIG. 1. The spin-up (top) and total ((a) in bottom) spec- 
tral function estimate for q^H/Tq = 48 {B/c = 1.1). For com- 
parison to the experimental curve (d) of Goldhaber-Gordon 
et al. [16] in a field of 7.5 T from Vds = ~2Q0^V to 200^iV 
with g^H/To « 50, (b) is (a) broadened by a temperature 
T = 90 mK [16] and (c) is the EOM result broadened by the 
same amount. 

The spectral function estimate for spin-up electrons 
follows from the single-spinon density of states by taking 
p-\{uj) = i/^crdwl) with a = +^ for w > and a — ~\ 
for w < 0, where a is the spin relative to the ground 
state. States of spin above the ground state can be 
obtained equivalently by adding a spin-up electron or by 
subtracting a spin-down electron, so the Kondo part of 



the spectral function is always symmetric under the com- 
bined transformation H — > — -ff, E —E. 

The shape of the spectral function is determined by 
the ratio gfiH/To- Fig. 1 shows the spin-up and total 
spectral functions for gfj,H/To = 48. The peak in the 
single-spin spectral function has width 2Tq at zero field 
and shifts away from the Fermi level and broadens with 
increasing field (Fig. 2). However, the width of the peak 
decreases as a fraction of gfiH. The shift of the peak is al- 
ways greater than the noninteracting level value g^H/2, 
starting at about two-thirds giiH for small gfj,H/To and 
rising slowly to gfiH as gfiH/To — > oo. The effect of 
temperature on the equilibrium spectral function will be 
slight until T > Tq; the effect of high temperature will 
be to blur of the Kondo correlations (and hence broaden 
and shift of the peaks toward the EOM result), but this is 
a smaller effect experimentally than the nonequilibrium 
broadening. The behavior of the finite-T magnetization 
suggests that the temperature needs to be a few times Tq 
for a strong effect. 





FIG. 2. The peak shift (top) in units of gfiH and broad- 
ening (bottom) in units of 2To, for various values of the ratio 
gfj,H/To. The dotted lines are limiting forms (12) and (13). 

The spectral function estimate can be compared to 
existing measurements of the H > differential con- 
ductance through a Kondo impurity, p^-p^ Within the 
NCA at = 0, the spectral function for nonzero bias 
Vds is essentially two broadened copies of the Vds ~ 
spectral function, one shifted up by Vds/2 and one down 
by Vds/2. ^ Note that Vds is the applied (drain-source) 
voltage across the dot, rather than the gate potential. 
The thermal broadening (T — 90 mK) and additional 
nonequilibrium broadening in [l^ mean that the ob- 
served peak is wider than the zero-temperature predic- 
tion, but the inward shift of the peak center should still 
be visible, and in and a smaller splitting than the 
EOM result of |^ was observed. Fig. 1 shows the agree- 
ment of our calculated peak location with that measured 
in Ina for \g\ = 0.36, the value suggested by electron 
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spin-resonance data on 2DEGs for the 7.5 Tesla field 
in that experiment. The value \g\ = 0.30 required for 
agreement with the EOM result is less likely. The ex- 
perimental peaks in Fig. 1 show some nonequilbrium 
broadening ||], but the parameters needed to estimate 
this broadening were not measured; convolving the theo- 
retical curve with a Lorentzian of variable width gives a 
very good fit to experiment. 

Because g in heterostructures is difficult to measure 
directly, a more convincing experimental demonstration 
of the splitting could be obtained by varying g^H/Tg. 
In experiment the nonequilibrium broadening is rel- 
atively large and may blur the correlations enough that 
the EOM result becomes apphcable. 

An asymptotic formula for the peak shape in the limit 
of large magnetic field can be derived from the largc-B 
limit of d) and CTgs(A) « 2c/'n:{AX^ + c^) : 



^ = l-^+0{x-Hog{x)), 
gjJLti Znx 



= B-X. (10) 



The result for the peak location in this limit is 



T" 



T" 



giiH 



2ttB 



2\og{g^iH^e/ToV2^) 



In the large- _ff limit the split between spin-up and spin- 
down peaks in the spectral function is 2gfiH since spin-up 
and spin-down peaks each shift by g^iH . The relative cor- 
rection in ([III ) is half that in the impurity magnetization 
M « {gtiH /2)[1 - l/\og{g^lH^e/ToV2^)]. 

The peak width (FWHM) in the hmit > To is 



g^iH 



27rB2 4 log2 {gfiH^/ToV2^) ' 



(12) 



The narrowness of the peak compared to giiH in the 
large- iJ limit is consistent with the EOM results. |^ 

There are two subtleties in the calculation worth men- 
tioning. The above expression for the magnetization 
gives a value continuously varying from AM = 1/2 for 
hole-like excitations with B to AAI = 1 for particle- 
like excitations with X <^ B. There are not continuous- 
spin excitations; rather the physical excitation combines 
a dressed spinon with an infinitesimal shift of the chem- 
ical potential for spinous B to give an integer or half- 
integer total spin. The same phenomenon of apparently 
continuous quantum numbers was seen in the charge sec- 
tor of the asymmetric Anderson model. At nonzero 
H and T = there is a small jump in the single-spinon 
density of states at E = giJ.H . This discontinuity is bal- 
anced by the appearance of new types of excitations, so 
that the total (thermodynamic) density of states is con- 
tinuous but for a (5-function at E = giJ-H from a state in 
the same SU (2) multiplet as the ground state. 

We now suggest a reason why large- iV techniques like 
the NCA give a spurious peak at the Fermi level for fi- 
nite H. Recall that the impurity orbitals transform under 



SU{N) rather than in a high-S representation of SU(2). 
The N degenerate orbitals split in a generalized mag- 
netic field as follows: A^/2 shift up by gfiH/2 and N/2 
shift down by the same amount. |^ The — > oo mag- 
netization curves for this and other splittings are derived 
in ||2^ . The Kondo peak near the origin then arises from 
transitions between states at the same energy, while the 
split peaks arise from transitions between states at dif- 
ferent energy. Counting the weight of each peak in this 
argument gives the correct value of the spin susceptibil- 
ity Wilson ratio, N/{N - 1). [§ The A^ = 2 problem 
is qualitatively different than the A^ > 4 problem since 
there are no degenerate orbitals after the magnetic field 
is applied, so it is not surprising that the analytic con- 
tinuation down to A^ = 2 gives a false peak at the origin, 
which is a real peak for A^ > 4. The problem of the con- 
ductance through a quantum dot in the Kondo regime is 
a fine example of the subtle properties of nonequlibrium 
transport in strongly correlated systems. 

J. E. M. acknowledges support from the Hertz Foun- 
dation and X.-G. W. from NSF grant DMR98-08941. 



Note added in proof: A preprint by T. Costi (cond 
mat/0004037[ ) finds comparable results for H <To from 



an NRG calculation of the spectral function. 
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